q-state Potts model on the Apollonian network 
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The q-state Potts model is studied on the ApoUonian network with Monte Carlo simulations 
and the Transfer Matrix method. The spontaneous magnetization, correlation length, entropy, and 
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scaling functions in temperature and q are proposed. A quantitative agreement is found between 
results from both methods. No critical behavior is observed in the thermodynamic limit for any 
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I. INTRODUCTION 



In the last five years, the Apollonian network (AN) 
[T] has attracted a lot of attention from the community 
working on models on non-Euclidean lattices. In partic- 
ular complex networks j2n4j found widespread use in the 
investigation of most diverse scientific topics, as they can 
be able to represent connections between individual de- 
grees of freedom in many complex systems. Also the be- 
havior of magnetic models on random complex networks 
has been investigated [3 [6]. AN's have the appealing 
advantages of geometrical sets defined through an exact 
inflation rule, where renormalization techniques, leading 
to exact results can be applied [ZHl]. In fact, this type 
of techniques has been considered for different hierarchi- 
cal structures to study critical phenomena [TDHIl]. This 
important property has also motivated using AN, as a 
first approximation, to study problems from several dif- 
ferent areas, such as physiology, geology, communication, 
energy, and fluid transport [T5H^ . 

The AN geometrical construction can be obtained re- 
cursively by initially taking three nodes in the vertices of 
an equilateral triangle, as shown in Fig. [ll A new node 
is then inserted in the center linked to those three. Se- 
quentially, new nodes are included, linked to each set of 
three connected nodes [TJ [5S] . The resulting network is 
scale free (power-law distribution of node degrees) and 
satisfies basic features of small- world networks, like large 
clustering coefficient and average minimal path i ~ InN, 
where N represents the number of nodes in the network. 

In this work we concentrate on the properties of ferro- 
magnetic Potts model on AN's. The Potts model [271 [^ 
represents a natural extension of the binary Ising model, 
where each spin variable is allowed to occupy a larger 
number q of independent states. This well known feature 
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FIG. 1: Apollonian network of generation two. 



is sufficient to produce a richer spectrum of phenomena 
as compared to the binary Ising model, as the influence 
of q on the properties of the phase transitions for models 
defined on Euclidean lattices. In turn, this raises the im- 
mediate question: how the Potts model, for g > 2, differs 
from the Ising model (g = 2) on AN's, where no phase 
transition at finite temperature has been detected. Does 
an increase in the value of q lead to a different scenario? 
The investigation of Ising models has already con- 
sidered a larger number of different situations. These 
include the next-neighbor interaction, which can be of 
ferro- or antiferromagnetic nature, coupling constants de- 
pending upon the generation where they were introduced 
into the model, or on the number of neighbors of each site 
p9l [30] . Also the effect of quenched disorder on the be- 
havior of Ising variables has been analyzed |15] . In this 
work, however, we consider only the simplest situation, 
i.e., homogeneous ferromagnetic couplings with a homo- 
geneous magnetic field. Our results are based on the 
independent use of two techniques: the standard Monte- 
Carlo (MC) simulations and the transfer matrix (TM) 
technique to numerically compute the partition function 
of the system. In both situations, we analyzed how the 
results depend on the size of the network. MC simu- 



lations consider up to 9844 lattice sites. On the other 
hand, TM technique aUows to go to larger sizes. In fact, 
the size can be chosen always sufficiently high to require 
a numerical convergence of the thermodynamical func- 
tions. As we show here, quite good agreement between 
the results from the two different techniques can been 
achieved. 

The paper is organized as follows: Section In] intro- 
duces the basic properties of AN networks and of the 
Potts model. Section [llT| brings some details of the used 
TM method. Section ITVl is divided into two subsections 
where we discuss, separately, the results obtained from 
MC simulations and TM maps. Finally, in Section \V\ we 
present our concluding remarks. 



II. APOLLONIAN NETWORK AND POTTS 
MODEL 



where the sum runs over directly connected pairs ij, 
yielding nearest-neighbor interactions, and the delta 
function, S (ct^, CTj), is unity when i and j are in the same 
state and zero otherwise. Jij is the coupling constant 
which, for simplicity, we consider to be the same for all 
interconnected pairs, Jij = J, and we take the limit 
J/f^B = 1, where kg is the Boltzmann constant. We 
use here the language of spins but, in fact, the q-state 
Potts model can also be applied to gauge theory, biolog- 
ical patterns, opinion dynamics, and image processing 
P71 [5T1 132] . More recently, a generalized version of the 
model has been proposed to study the topology of net- 
works through the identification of different communities 
as well as the overlap between them [331 - 135] . In the next 
section we discuss how to use a transfer matrix formal- 
ism to study the q-state Potts model on the Apollonian 
network. 



The simplest Apollonian network (AN) [2 [5S] is ob- 
tained recursively by first placing three nodes at the cor- 
ner of a triangle (generation 0) . A new site is put into the 
triangle and connected to all three corner nodes forming 
three new triangles (generation 1). Then, at each gener- 
ation, a node is placed into each triangle and connected 
with its three corner nodes. Being n the generation, the 
number of nodes, N, is given by 



N{n) = 



3"-f5 



and the number of edges, E, by 



E{n) 



gn+l 



(1) 



(2) 



In the limit of large n, it is straightforward that 
E{n)/N{n) — >■ 3, i.e., on average each node is linked 
with six other nodes. The distribution of links is het- 
erogeneous and the network is scale free, with a degree 
exponent 7 « 1.585 [T. 

In the q-state Potts model, each node of the network 
contains a spin which can assume q different states, cr, 
i.e., a = l,2,3,...,q. The Hamiltonian of the model is 
then. 



"H = - ^ Jij(5 {a„ aj)- h^cTi 



(3) 



III. TRANSFER MATRIX AND RECURRENCE 
MAPS 

We have used a transfer matrix (TM) formalism to nu- 
merically evaluate the partition function for several Ising 
models on Apollonian networks [29l [30] . This is a very 
useful method as it yields the properties of the system 
for any given generation n. It is also possible to reach 
numerically the thermodynamic limit, where the free en- 
ergy per spin and its derivatives can be obtained to any 
pre-established precision (usually < 10~^^). The method 
takes advantage from the AN scale invariance when we 
go from generation n to n -I- 1, and from the fact that 
partial sums over all spins at generation n can be recur- 
sively performed when we write the partition function for 
generation n + 1. 

Of course the same strategy can also be used when 
Potts spins are considered, the only difference being that 
we must consider matrices of order q x q and q x q^. For 
the sake of simplicity, let us write down explicitly the 
case g = 3, with h ~ 0. The generation zero, n — 0, 
consists only of the three spins placed at the vertices of 
the largest triangle in Fig. [Tl If we perform a partial 
trace over the spin on the lower vertex, the interaction 
between the sites i and k can be condensed in a single 
TM as: 



'ao bo bo 
Mo = I 60 flo bo 
,bo bo oq. 



b{2ab 
\b{2ab 



(q- 


-l)fe2) b{2ab + {q - 2)b^) b{2ab+{q- 


-2)62 


{q- 


-2)fe2) a{a^ + {q-l)9) b{2ab + [q - 


-2)6" 


(q- 


-2)fe2) b{2ab+{q-2)9) aia^ + {q- 


- ^W] 



(4) 



r 



where a = exp(/3J) and 6 = 1. To proceed further with 
the method and consider generation n — 1, it is necessary 



to define a.qxq^ TM Lq, which describes the interactions 
among sites i, j, and k. We use a column label k that 



depends on the pair (j, k) according to the lexicographic 
order, i.e., k = q{j — 1) + fc, so that 



^■o So So So So to So to So\ 

-^0 = I So So to So To So to So So = 

So to So to So So So So roj 

a3 ab^ 63 ab^ a9 



a^ ab"^ ab^ ab"^ 
ab"" ab"" 63 a62 
,a62 63 a9 b^ 



(5) 



ab"^ ab'^ ab"^ ab"^ 



To evaluate the magnetic properties, it is necessary 
to break the symmetry among the q states and insert 
a non zero field h ^ along one of the a (say a = 1) 
directions. This changes the simple structure of matrices 
M„ and L„, that then have a much larger number of 
matrix elements. By performing an explicit calculation 
of the matrix elements for larger values of q, it is possible 
to complete the set of maps used in this work, as listed 
in the Appendix. 



As discussed in detail in Rcf. 1 , transfer matrices Mi 
and Li can be expressed in terms of Mq and Lq as 



9 <7 






and 



9 



{Ll)ijk — 7 ,(-^o)ij£(-^o)t/fc(-^o)fcj£ 



(7) 



Note that the above matrix maps require that the ma- 
trix element to (t„) needs to be introduced only for in- 
teger q > 2. The matrix elements a„ and 6„ can be 
expressed as a„ — rn + {q—l)sn and 6„ = 2s„-|-((7 — 2)t„. 

Since the network grows according to a generation in- 
dependent inflation rule, Eqs. m and ([7| apply for any 
value of n, just by replacing 1 by n -I- 1 and by n. Such 
general matrix maps can be rewritten in terms of maps 
for the matrix elements of Ln+i in terms of those of L„ 
as 



r„+i = rl + {q- l)sl 

Sn+i = r„sl + sl + {q- 2)sntl 

t„+i = 3sltn + iq- 3)t3 



(8) 



from which the elements a„+i and 6„+i can be obtained. 
It is possible to directly evaluate the free energy and 
all other thermodynamic functions in terms of a„ and 
6„ or, equivalently, from r„, s„, and t„. Due to the 
M„'s particular form, its eigenvalues can be easily eval- 
uated as A„ = an + {q — 1)6„ and the (g^l)-degenerated 



Xn 



- bn- Since the numerical values of a^, 6„, r„, Sr, 



and t„ grow exponentially, it is convenient to write down 
recurrence maps that avoid numerical overflows when 
they are iterated. Since r„ is the fastest exponentially 
growing variable, this can be accomplished by deriving 
maps for the free energy /„ = ~T ln{An) / N {g) , the cor- 
relation length ^n = 1/ ln(A„/A„), and the auxiliary vari- 
able yn = tn/rn- An alternative definition would be to 
replace the map for ^„ by that for Xn = Sn/fn, from 
which the value of ^„ can be evaluated. This set of maps 
can be enlarged by working out explicit recurrence rela- 
tions for the derivatives of fn{T) and the other variables 
[yn, Xn, or £^n) with respect to the temperature. This 
way, the explicit temperature dependence of the entropy 
s{T) and the specific heat c{T) can be obtained. 



IV. RESULTS AND DISCUSSION 

To study the q-state Potts model on the Apollonian 
network we consider two different approaches: Monte 
Carlo simulations (MC) and the Transfer Matrix method 
(TM). Monte Carlo simulations are limited to small sys- 
tem sizes, so results are affected by finite-size effects. On 
the other hand, with the Transfer Matrix formalism large 
system sizes can be considered allowing to compute the 
thermodynamic and magnetic properties, numerically, in 
the thermodynamic limit. This section is then divided 
in two subsections. In the first we discuss the results 
obtained with Monte Carlo simulations for the magneti- 
zation, r7i(g,T), and the specific heat, c{q,T), as a func- 
tion of the temperature, T, and the number of states, 
q. A comparison between MC and TM results is also 
included. In the second part, not only the magnetiza- 
tion and specific heat but also the entropy, s{q,T), and 
the correlation length, ^(g, T), are obtained with the TM 
technique for larger system sizes. The thermodynamic 
limit is then discussed. 



A. Monte Carlo simulations 

We carried out Monte Carlo simulations for two differ- 
ent generations of the Apollonian network, 7 and 9, corre- 
sponding, respectively, to 1096 and 9844 spins (Eq. Q). 
For each system size, we considered two values of q (num- 
ber of states), 2 and 3, and we computed the magnetiza- 
tion and specific heat as a function of the temperature. 
For the Potts model jM] , the magnetization is defined as 



m{q,T) 



q{ni{q,T) - 1) 
q-1 



(9) 



where ni is the fraction of spins in the state 1, and the 
specific heat is defined as 



ciq,T) 



1 



J 



N V /crT 



[<u^ >~ <u>^] , (10) 



where N is the total number of spins and < u > and 
< u^ > are, respectively, the first and second moment 
of the energy per spin. All Monte Carlo results were 
averaged over 10^ samples. 

In Fig. ^ we show the magnetization, 'm{q,T), as a 
function of the temperature for q = 2 and q = 8, for 
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FIG. 2: (Color online) Magnetization, m{q,T), as a function 
of the temperature, T, for generation 7 (inset) and 9 (main 
plot). The considered number of states, q, are 2 (red squares) 
and 8 (blue circles). The lines were obtained with Transfer 
Matrices and the points with Monte Carlo simulations. Sim- 
ulational results were averaged over 10* samples. Error bars 
are within the point size. 
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FIG. 4: (Color online) a) Magnetization, m{q, T) as a func- 
tion of the temperature, T, for q = {2, 4, 8, 16, 32} (from right 
to left). The magnetization vanishes exponentially with the 
temperature (see inset) and the argument of the exponential 
decreases linearly with q. b) Correlation length, ^{q,T), as 
a function of the temperature, T, for the same values of q. 
Results obtained for generation 100. 



FIG. 3: (Color online) Specific heat, c{q,T), as a function of 
the temperature, T, for generation 7 (inset) and 9 (main plot). 
The considered number of states, q, are 2 (red squares) and 8 
(blue circles). The lines were obtained with Transfer Matrices 
and the points with Monte Carlo simulations. Simulational 
results were averaged over 10* samples. Error bars are within 
the point size. 



generation 9. The inset shows the same functions for 
generation 7. The points correspond to Monte Carlo re- 
sults and the lines to numerical results obtained with the 
Transfer Matrix method. An exponential decrease of the 
magnetization is observed. The larger the number of pos- 
sible states the steeper the change with temperature. We 
find an agreement between the Monte Carlo results and 
the ones from the Transfer Matrix method. 

Figure [3] shows the specific heat, c{q,T), as a function 
of temperature. Like in Fig. [2) two different values of 



states (2 and 8) and generations (7 and 9) are consid- 
ered. For both values of g a Schottky maximum in the 
specific heat is found for values of T below 2. With in- 
creasing number of states the maximum becomes sharper 
and occurs at lower temperature (discussed in the next 
subsection). Results for different generations reveal no 
significant size effects in the specific heat. A quantita- 
tive agreement between the Monte Carlo (points) and 
Transfer Matrices (lines) results is obtained. 



Results from Monte Carlo simulations are affected by 
finite-size effects. In the next subsection, we consider the 
Transfer Matrix formalism to study larger system sizes 
and discuss the behavior of the system in the thermody- 
namic limit. 
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FIG. 5: (Color online) Entropy, s{q,T), as a function of the 
temperature, T, for q = {2, 4, 8, 16, 32} (from bottom to top). 
Results obtained for generation 100. For all values of q the 
entropy converges to a fixed value for high temperatures. In 
the inset we see this value as a function of q at two different 
temperatures: 100 (red squares) and 1000 (blue dots). Above 
a certain temperature, the entropy scales logarithmically with 

q- 



B. Transfer Matrix 

We evaluated the partition function for the q-state 
Potts model on an Apollonian network by numerically 
iterating the transfer matrix maps introduced in Sec. |III| 
and detailed in the Appendix. We study thermodynamic 
and magnetic properties like the spontaneous magneti- 
zation, TO, the specific heat, c, the entropy, s, and the 
correlation length, ^, for different values of q, as a func- 
tion of the temperature, T. 

Figure WTa) shows the magnetization, rn{q,T), as a 
function of the temperature, T, for different values of the 
number of states, q = {2,4,8,16,32}. A smooth decay 
of the magnetization with the temperature is obtained, 

m(r,g)^exp(~0(g)T) , (11) 

where (j){q) is a linear function of q, 

(j){q) ^aq + b , (12) 

with units of the inverse of temperature. We have es- 
timated a = 0.18 ± 0.01 and b = 0.17 ± 0.03, and no 
significant finite-size effects are observed. For q — 2, 
the q-state Potts model is equivalent to the Ising model 
and we recover the previously reported behavior for this 
system [33]. Note that, the temperature in the Ising 
model, Tjsing, differs from the one of the Potts model, 
Tpotts, by a factor of 1/2, i.e., Tjsing = Tpottsf^- Conse- 
quently, a decay of the magnetization, in the Ising model, 
as exp(— T), corresponds to exp(— T/2) here. 

The behavior of the magnetization shows no evidence 
of a sharp order-disorder transition, which is not the case 



for the correlation length, $,{q,T), in Fig. |4Fb). For all 
values of q there is a well-defined temperature, T*, below 
which the correlation length numerically diverges when 
compared with the value of ^ in the disordered phase. 
Here we chose the divergence threshold to be 10'^®, but 
any other value would lead to the same qualitative pic- 
ture. Above T* the correlation length attains finite val- 
ues. This apparent transition is in fact a finite-size effect. 
As observed for the Ising model [25] , the value of T* goes 
linearly with the generation, n, being infinite in the ther- 
modynamic limit, i.e., T* ^ 9{q)n, where, for large q, 



9{q) ^ log" (9) 



(13) 



and T] — —3.81 ± 0.07. The value of T* decreases with q 
according to. 



T*r^logHq) 



(14) 



with ^ — —2.74 ± 0.04. The absence of criticality is in 
agreement with previous results, where several models 
characterized by criticality on periodic networks show 
no critical transition, for finite temperature, in scale-free 
networks with a degree exponent 7 < 3 [TTHID] . In this 
regime, an ordered phase is observed at any tempera- 
ture. Kaplan, Hinczewski, and Berker [22] reported that 
in the presence of quenched disorder the ordered phases 
are still robust and persist for the entire range of disor- 
der. Recently, Igloi and Turban [37 have considered a 
mean-field version of the Potts model and reported that 
an order/disorder transition solely occurs for 7 > 3. A 
(7(7) can then be defined above which the order of the 
transition changes from second- to first order, like on pe- 
riodic lattices. 

In Fig.lslwe draw the curves of the entropy, s{q, T), as a 
function of the temperature. For all considered values of 
q, the entropy rises with T and saturates at a maximum 
for large temperatures (above 10). The inset of Fig. [s] 
contains the entropy as a function of the number of states, 
q, for two different temperatures, 100 and 1000. For large 
T, the entropy depends logarithmically on q, 



siq,T)^log{q) 



(15) 



These results have no dependence on the system size. As 
for the specific heat (discussed below), the curve of the 
entropy converges fast with the generation number. 

The specific heat, c{q,T), is shown in Fig. [6]^ a). As al- 
ready observed for smaller system sizes with Monte Carlo 
simulations, the dependence of the specific heat on the 
temperature is smooth with a peak at a temperature be- 
low 2. No divergence of the maximum with the system 
size is observed as would be expected for a transition. 
In fact, this peak corresponds to a Schottky maximum 
and is due to a fast increase in the entropy (see Fig. I5| . 
The larger the value of q the sharper the peak and the 
lower the temperature at which it occurs. For large T the 
specific heat diminishes with temperature according to a 
power law with an exponent 2.00 ± 0.02. This exponent 
is the same for different values of q and the behavior is 
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FIG. 6: (Color online) (a) Specific heat, c(q, T), as a function 
of the temperature, T, for q = {2, 4, 8, 16, 32} (from bottom 
to top). Results are obtained for generation 100. For high 
temperatures, the specific heat scales with the temperature 
according to a power law (inset in (a)), (b) Specific heat as 
a function of q at two temperatures: 100 (red squares) and 
1000 (blue dots). 



dynamic limit for any value of the number of states, q. 
The specific heat is characterized by a Schottky maxi- 
mum which becomes sharper with increasing q, without 
divergence as expected for a transition. As previously re- 
ported for the Ising model on the same network [^ [SU] , 
the specific heat converges rapidly with the generation 
number to the thermodynamic limit. 

In the present work we consider spins interacting ferro- 
magnetically with its nearest neighbors in the absence of 
a magnetic field. A more general version of the model can 
be considered where pairs of spins interact with further 
neighbors, e.g., next-nearest neighbors, and a magnetic 
field, either uniform or site dependent. Besides, it would 
be interesting to analyze the effect of different coupling 
constants between spins based on their generation [30] . 
to observe possible occurrence of a critical behavior in 
the thermodynamic and magnetic properties when the 
temperature changes. The considered methodology can 
also be taken to study the properties of the recently pro- 
posed versions of the model to identify communities in 
networks I33H351. 



VI. APPENDIX 

The map for the free energy assumes a simpler 
form if we consider an alternative definition / — 
— Tln(r„)/7V(g), from which the value / can be easily 
obtained. In fact, in the limit of large n, f and / become 
identical. The full set of maps, which depends also on 
the variable g — exp{l3h) with /? = 1/r, reads: 



/n+i = ^^ - T^ H9 + {<!- l)vl) (17) 



Nn+l Nn+l 



^'n+i ul + gwl + {q- 2)xl 

Un+l — — —, TT-^ (18) 



rn+1 



g + {q- l)w3 



independent of the system size. In this power-law regime, 
for a fixed temperature, the specific heat scales linearly 
with q (FigjoFb)). Therefore, for large T, the following 
scaling law can be postulated, 



c{q,T)^l/qT- 



(16) 



s„+i gvl -f v^wl + {q- 2)vnyl /, n\ 

"^"+1 = I — = rm — t::3. (19) 



rn+l 



g+{q~ l)v^ 



_ Sn+i _ gvlwn + Unwl + {q - 2)xnyl /on^ 

■Tn+l 9+{q- 1)W„ 



V. CONCLUSIONS 

We studied the q-state Potts model on the Apollonian 
network through Monte Carlo simulations and Trans- 
fer Matrix method. Different scaling relations were ob- 
tained for magnetic and thermodynamic properties like 
the spontaneous magnetization, correlation length, en- 
tropy, and specific heat. We have shown that the magne- 
tization decays smoothly with the temperature and the 
greater the number of states the steeper the decay. How- 
ever, no order-disorder transition is found in the thernio- 



^n+l 



Sn+1 _ gwnvl + ^nxi + x^j + (^ - 3)a;„z2 



rn+l 



g + {q- l)vl 



n+l 



+ {q- l)w3 



?'«+! 



g + {q- l)v^ 



(21) 



tn+l gVnvl -f iWnXnVn + [q ^ 3)2/,^Z„ 

Vn+1 = z = r-TZ TTZ^ (22) 



_ <«+i _ gyl + 3xlz„ + {q- 4)4 

^n+l — z ^ _ , /^ TVZ^ v^'^/ 



New variables r„,s„,Sn, and !„ have been introduced 
because the external field h reduces the problem sym- 
metry, as reflected in a larger number of distinct matrix 
elements. Note that fi = ri,si = si = si, and ti = ti. 
If /i = 0, the number of independent equations is reduced 
to three. 
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